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I.  INTRODUCTION 


Adaptive  arrays  based  on  the  LMS  algorithm  of  Widrow,  et  a1[l]  can  be 
used  to  protect  communication  systems  from  interference.  These  antennas 
optimize  reception  of  a  desired  signal  while  simultaneously  nulling 
interference.  For  a  given  set  of  incoming  signals,  the  LMS  array  yields  the 
maximum  attainable  signal -to-interference-plus-noise  ratio  (SINR)  at  the  array 
output[2,3]. 

Most  studies  of  adaptive  array  performance  have  assumed  interference 
signals  that  have  constant  power.  However,  a  modulated  interference  signal 
may  be  more  difficult  for  the  array  to  null.  If  the  modulation  rate  of  the 
interference  is  close  to  the  natural  response  rate  of  the  array,  such  a  signal 
may  keep  the  array  in  a  continual  transient  state. 

In  a  previous  paper[4],  one  of  the  authors  examined  the  performance  of  an 
LMS  array  with  pulsed  interference.  It  was  shown,  for  example,  that  when  the 
array  is  used  to  receive  a  differential  phase-shift  keyed  (DPSK)  signal, 
pulsed  interference  can  increase  the  bit  error  probability  more  than  CW 
interference  does.  However,  the  increase  due  to  the  pulse  modulation  is 
small,  and  occurs  only  if  the  pulse  repetition  rate,  duty  cycle  and 
interference  power  are  all  properly  chosen.  Otherwise,  pulsed  interference 
has  less  effect  than  CW  interference. 

In  this  report,  we  study  the  effect  of  another  type  of  envelope  modulated 
interference  on  the  LMS  array.  Specifically,  we  consider  an  interference 
signal  with  double-sideband,  suppressed  carrier,  sinusoidal  envelope 
modulation.  This  type  of  signal  has  been  chosen  for  two  reasons.  First, 
such  interference  contains  a  spectral  line  on  each  side  of  the  desired  signal 
frequency;  by  changing  the  interference  modulation  frequency,  we  can  explore 
the  frequency  response  of  the  LMS  loop.  (The  LMS  loop  is  nonlinear  in  the 
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input  signals.)  Second,  such  modulation  yields  a  system  of  differential 
equations  for  the  array  weights  that  can  be  solved.  In  general,  an 
interference  signal  with  arbitrary  modulation  leads  to  an  intractable 
mathematical  problem. 

As  we  shall  show,  modulated  interference  of  this  type  has  two  effects  on 
the  array.  First,  it  causes  the  array  to  modulate  the  desired  signal  envelope 
(but  not  its  phase!).  Second,  it  makes  the  output  S I NR  from  the  array  vary 
with  time.  However,  it  turns  out  that  this  S I  NR  variation  causes  almost  no 
change  in  average  bit  error  probability  (as  compared  with  CW  interference) 
when  the  array  is  used  in  a  digital  communication  system. 

In  section  II  of  the  report  we  establish  notation,  define  the  interference 
signal  and  solve  for  the  weight  response  of  the  adaptive  array.  In  section 
III  we  present  calculated  results  obtained  from  the  equations  in  section  II 
and  discuss  the  effect  of  each  interference  signal  parameter  on  the  array 
performance.  Section  IV  contains  the  conclusions. 

II.  FORMULATION  OF  THE  PROBLEM 

Consider  an  LMS  adaptive  array[l|  consisting  of  three  isotropic  elements 
a  half  wavelength  apart,  as  shown  in  Figure  1,  The  analytic  signal  xj(t)  from 
element  j  is  multiplied  by  complex  weight  wj  and  then  summed  to  produce  the 
array  output  s(t).  The  error  signal  c(t),  which  is  the  difference  between  the 
array  output  signal  and  the  reference  signal  r(t),  forms  the  input  to  a 
feedback  system  that  adjusts  the  wj .  In  the  LMS  array,  the  weights  satisfy 
the  system  of  equations[l ,5] 

$£•  +  k*W  =  kS,  (1) 

where  W=(w^,  W2,  wg)T  is  the  weight  vector,  <$  is  the  covariance  matrix. 


?. 


*  =  E[X*Xt], 


(2) 


S  is  the  reference  correlation  vector 


S  =  E[X*r(t)],  (3) 

and  k  is  the  LMS  loop  gain.  In  these  equations,  X  is  the  signal  vector, 

X  =  [x! (t ) ,  x2(t),  x3(t)]T,  (4) 

T  denotes  transpose,  *  complex  conjugate,  and  E[»]  expectation. 

We  assume  that  a  desired  and  an  interference  signal  are  incident  on  the 
array  and  that  thermal  noise  is  also  present  in  each  element  signal.  The 
signal  vector  then  contains  three  terms. 


X  =  Xd  +  Xi  +  Xn, 


(5) 


where  Xj,  Xj  and  Xn  are  the  desired,  interference  and  thermal  noise  vectors, 
respectively. 

Specifically,  we  assume  a  CW  desired  signal  is  incident  from  angle  9,j 
relative  to  broadside,  (e  is  defined  in  Figure  1.)  The  desired  signal  vector 
is  then 


v  a  (aot+’l'd)I, 

Xd  =  Ade  IJd, 


(6) 


where  Ad  is  the  amplitude,  <lq  is  the  carrier  frequency,  ipd  is  the  carrier 
phase  angle,  and  Ud  is  a  vector  containing  the  interelement  phase  shifts, 

Ud  «  (1,  e“j*d,  e-j2*<l)I  (7) 


with 

4>d  =  irsinOd.  (8) 

We  assume  ^  to  be  a  random  variable  uniformly  distributed  on  (0,?tt). 

Next,  we  assume  an  envelope  modulated  interference  signal  as  shown  in 
Figure  ?.,  arriving  from  angle  .  The  interference  signal  vector  is  then 


INTERFERENCE 

SIGNAL 


Figure  2.  Sinusoidally  Modulated  Interference 
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(9) 


ai  (t) 

x.  =  eJ(v>t+^)  ai  {t-T-|  )e"'ia®T1 
ai(t-2Ti)e"'i2‘,«Tl 

where  &]{t)  is  the  envelope  modulation  as  received  on  element  1,  ^  is  the 
carrier  phase  angle,  and  T-j  is  the  interelement  time  delay  given  by 

T-j  =  —  sinei .  (10) 

tK> 

We  assume  ^  is  a  random  variable  uniformly  distributed  on  (0»2ir)  and 
statistically  independent  of 

In  this  paper,  we  shall  study  the  case  where  a-j(t)  is  a  simple  sinusoid, 

a^  (t)  =  A-j cos^t ,  (11) 

with  Aj  and  u^,  the  amplitude  and  frequency,  respectively.  I.e.,  the 
interference  is  a  double-sideband,  suppressed  carrier  signal.  As  mentioned 
above,  we  have  chosen  this  type  of  modulation  for  two  reasons.  First,  this 
signal  has  two  spectral  lines,  one  on  each  side  of  the  desired  signal 
frequency.  By  varying  we  can  change  the  separation  between  the  desired  and 
interference  spectral  lines  and  examine  the  frequency  response  of  the  array 
weights.  Second,  this  type  of  modulation  leads  to  a  tractable  mathematical 
problem.  More  complicated  modulations  yield  differential  equations  for  the 
weight  vector  that  are  more  difficult  to  solve.* 

We  shall  assume  the  frequency  u^,  in  (11)  is  very  small  compared  to  the 
carrier  frequency  (lq.  Under  this  condition  the  time  delay  T-j  represents  a 
negligible  delay  in  the  modulation  envelope,  so  the  modulation  on  each  element 
signal  is  essentially  the  same, 

ai(t)  =  ai(t-Ti)  =  ai(t-2Ti).  (1?) 

*Specifical ly,  the  modulation  in  (11)  yields  a  3-term  recursion  relation  in 
Eg. (57)  below.  If,  for  example,  a-j(t)  is  assumed  to  contain  a  carrier  term  in 
addition  to  the  two  sidebands,  a  5-term  recursion  relation  results  instead  of 
(57).  The  weight  differential  equation  then  cannot  be  solved  by  the  methods 
presented  here. 


In  this  case  reduces  to 


Xi  =  Aicosuffltej^t  +  'i'^!Ji ,  (13) 

where 

Ui  »  [1,  e"j*,  e“j2*1]T,  (14) 

with 


<()-j  =  usinQ-j .  (IS) 

Finally,  we  assume  the  thermal  noise  vector  is  given  by 

Xn  =  [ni(t),  n2(t),  n3(t ) ]T ,  (16) 

where  the  nj(t)  are  zero-mean,  gaussian  thermal  noise  voltages,  all 
statistically  independent  of  each  other,  and  each  of  power  a2.  Thus, 


ECn j  (t )n)< (t ) ]~o26jk ,  (17) 

with  <5ij  the  Kronecker  delta.  Also,  we  assume  the  fij(t)  are  statistically 
independent  of  ^d  and  ip-j  . 

Under  these  assumptions,  the  covariance  matrix  in  (2)  becomes 

$  =  <J>d  +  (18) 

where 

*d  =  Ad2Ud*Ud',  (19) 

4>j  =  Aj  ^cos^omtlJ,  *U-j T 

=  (l/2)Aj2[l+cos?u111t]l)i*UiT,  (70) 

and 

*n  =  o2I.  (21) 


where  I  is  the  identity  matrix. 
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To  compute  the  reference  correlation  vector  S  in  (3),  we  must  first  define 
the  reference  signal  r(t).  In  practice,  the  reference  signal  is  usually  derived 
from  the  array  output  [6, 7].  It  must  be  a  signal  correlated  with  the  desired 
signal  and  uncorrelated  with  the  interference.  Here  we  assume  the  reference 
signal  to  be 

r  (t)  =  Are^<i°t+^c*^ .  {?.?) 

Eq.(3)  then  yields 

S  =  ArAdUd*.  (23) 

Eqs.(l8)-(21)  and  (23)  can  now  be  inserted  in  (1)  to  give  the  differential 
equation  for  the  weight  vector  W, 

^  +  k[a2I+Ad2Ud*UdT+(l/2)Ai2(l+cos2ufflt)Ui*UiT]W  =  kArAdUd*.  (24) 

Before  attempting  to  solve  (24),  it  is  helpful  to  put  it  in  normalized  form. 
First,  dividing  by  ko2  gives 

+  [I  +  CdUd*UdT+(l/2)^i(l+cos2afl,,t')IJi*UiT]W(f)  =  Ud*.  (25) 

where 

£d=  =  input  signal-to-noise  ratio  (SNR)  per  element, 

a2 

A-  2 

r-j  =  — l  =  peak  input  interference-to-noise  ratio  (INR)  per  element. 
a2 

and  where  we  have  also  used  normalized  time  and  frequency  variables, 
on,'3  ttm/k a2  =  normalized  modulation  frequency. 


t'=  ko2t  =  normalized  time. 


Next  we  note  that  the  constant  Ar/o  on  the  right  side  of  (25)  will  just  appear 
as  a  scale  factor  in  the  solution  for  W.  It  has  no  effect  on  the  array  output 
signal -to-noise  ratios  to  be  computed  below.  Hence  we  arbitrarily  set  Ar/o=l  to 
eliminate  it.  Finally,  by  defining 

*o  =  1  +  SdUd^i1  +  (l/2niUi*U1T,  (26) 

we  can  write  (25)  in  the  form, 

C4o+(l/2) CiCOS (2tim't '  )U-j*U-jT]W(t 1 )  =  /pd\  (27) 

Eq.(27)  is  a  linear  vector  differential  equation  with  a  constant  source 
term  on  the  right  but  with  periodic  coefficients.  The  general  properties  of 
such  differential  equations  have  been  described  by  O' Angel o[8].  He  has  shown 
that  the  solution  for  W(t')  will  be  a  periodic  function  of  time  once  the  initial 
transients  have  died  out. 

In  this  paper,  we  shall  not  consider  the  initial  transient  behavior,  but 
shall  concentrate  on  the  periodic  steady-state  solution  for  W(t').  Once  W(t') 
is  periodic,  it  can  be  represented  as  a  Fourier  series 

W(t')  =  l  Cnej2nafl’,t',  (28) 

n=-<» 

where  Cn  is  a  vector  Fourier  coefficient.  Substituting  (28)  into  (27), 
replacing  cos^u^'t')  wit*1  exponentials,  and  collecting  terms  with  the  same 
frequency,  we  find  that  the  coefficients  Cn  must  satisfy 

J.  C(«o+j2nofl,,I)Cn  +  (l/4)CiUi%-T(Cn.1+Cn+1)]ej2n^,t  =  /ejud*.  (29) 
n=-« 

Enforcing  this  equation  for  each  frequency  component  separately  gives,  for  n=0. 


«oco+(l/4)^i‘Ji*UiT(c.1+C1)  =  /^Ud*. 


(30) 


and  for  n*0. 


( <mti '  I Kn  *  ( Uj^(Cn_x+Cn-i-x)  =  0.  (31) 

These  equations  may  be  solved  for  the  Cn  by  expressing  each  Cn  in  terms  of 
its  components  parallel  and  perpendicular  to  the  vector  U-j*.  To  do  this,  we 
first  form  a  set  of  three  orthonormal  basis  vectors*  as  follows.  Let  ex  be  a 
unit  vector  parallel  to  U-j*, 

ex  =  ,  (32) 

Vuj* 

and  let  e2  and  e3  be  unit  vectors  perpendicular  to  each  other  and  to  ex,  i.e., 
ej*ek  =  5jk»  (33) 

(where  t  is  conjugate  transpose).  In  particular,  let  e?  be  in  the  plane  defined 
by  Ud*  and  Uj*. 

^2  =  c[Ud*-<Uj*],  (34) 

where  5  and  *  are  constants.  Enforcing  the  orthogonality  condition  ex^e2=0 
and  choosing  c  so  e2  has  unit  magnitude  yields, 

[Ud*  -  ^Mli*].  (35) 

UiTUi* 

63  is  easily  found  from  ex  and  e2  but  will  not  be  needed  in  the  sequel,  so  we 
shall  not  compute  it  explicitly. 

Each  coefficient  Cn  may  be  written  in  terms  of  the  unit  vectors  e^  as 
3 

^n  =  1  an,kek»  (36) 

k=l 

where  the  0^  ^  are  scalar  coefficients,  is  the  component  of  Cn  along  the 

unit  vector  e^. 

*The  array  has  3  elements,  so  W(t')  and  Cn  each  have  3  components.  Hence  3 
basis  vectors  are  needed  to  span  the  space  for  Cn. 


e2  = 


✓  UdTUd*- 

UiTUl 
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Substituting  Eq.(36)  into  (30)  and  (31)  and  utilizing  the  fact  that 

U-jTej=0  for  j*l,  we  obtain 
3 

4o  I  «0,kek  +  ( l/4)Ci  ttJi *1  (<M,l+al,l)  =  >^dud*»  (37) 

k=l 

and 

3 

l  on.k^o+j20^' 1  )ek+(  1/4) Ci  I U-j *  1  (on_i  x+an+i  i)=0,  n*0.  (38) 

k=l 

Multiplying  these  two  equations  by  ej'*'  on  the  left  gives 
3 

l  oO.kSjk  +  n/4)^  I  Ui  *  |  (a_ifi+aifi)6ji  =  /|dejtud*»  (39) 

k=l 

and 

3  *  ? 

I  an>k(ejk+j2naffl,«jk)  +  (l/4)?i|Ui  |  (an-l,l+o«n+1>1)6j1=0,  n*0,  (40) 

k=l 

where  we  have  defined 

Sjk  =  ejfMk-  (41) 

Eqs . (39)  and  (40)  hold  for  l<j<3.  THe  values  of  8jk  are  readily  found  by 
substituting  (26),  (32)  and  (35)  into  (41)  and  making  use  of  (33)  for  eg. 
results  are 


Bn  =  1+  2^(  1  /2 )  Ui  TUi  * , 

UiTUi* 

(42) 

822  =  l+CdP» 

(43) 

Si?  =  8?i  -  ^  -1"  d~. 

^Ui TUi * 

(44) 

Bl3  =  S31  =  823  =  S32  =  0, 

(45) 

833  = 

(46) 

The 
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where 


(47) 


In  addition,  one  finds 


«lV  - 

(48) 

V^i* 

e2tUd*  =  /p. 

(49) 

e3%*  =  0- 

(50) 

Our  problem  now  is  to  solve  (39)  and  (40)  for  the  We  may  do  this  as 

follows.  Writing  out  (39)  and  (40)  for  l<jO  gives  the  six  equations, 


3liaO,l  +  012oO,2+(1/4)^i  |Ui*|  2(a-14+<n,l)  =  /Sdelt|Jd*»  (51) 

®21aO,l+®22aO,2=,^de2^^d*»  (52) 

a0,3=0>  (53) 

(611+320^'  )an>1+Bi2an>2=-(l/4)ri  | Ui*|2(cjn_i>i+an+itl),  (54) 

821an,l+(fl22+j2n(Mn'  )°>n,2=0»  (55) 

and 

<*n,3  =  °*  (56) 

where  the  last  three  equations  hold  for  n*0.  Solving  (54)  and  (55)  for  ant\  and 
«n,2  yields 

<*0,1  =  -Yn(an-1  ,l+an+l,l)  (57) 
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where 


(B??+j2nam‘ )  lUj*| 2 _ _ 

4[(0ii+j2n<Mnl  )(022+j2n(Hn')-Ul2|2l 


(59) 


Eq.(58)  allows  one  to  calculate  0^2  from  cy,^,  so  the  problem  is  reduced 
to  determining  the  To  find  the  cy,^,  we  must  solve  (57). 

Eq.(57)  is  a  three-term  recursion  relation  for  the  coefficients  an>i. 

Given  cy^i  for  two  successive  n,  one  could  calculate  all  the  remaining  an>i  from 
(57).  However,  if  one  starts  from  arbitrary  initial  values  for  the  an>i,  one 
finds  that  in  general  oy,^  will  not  approach  zero  as  n+t®,  so  the  series  in  (28) 
will  not  converge.  Only  the  correct  initial  values  of  «n>i  will  cause  the 
series  to  converge.* 

To  determine  the  correct  initial  values  for  the  cy,^,  we  manipulate  (57) 
into  two  forms. 


and 


-  «hJL . 
on-l.l 


-1 

IV  -9n.il  a 
7n  <*1,1 


for  n>0. 


(60) 


-A»l. 
an+l  ,1 


■j - - -  for  n<0. 

I  +  «n-l.l 
Tn  afl,l 


(61) 


Consider  first  the  behavior  of  «n>i  for  large  positive  n.  Eq.(59)  shows  that 
for  large  n. 


I  *  ,  MUianl  . 

Yn  ^ lUi*|2 


(62) 


*The  situation  here  is  similar  to  that  encountered  in  solving  the  Mathieu 
equation[9]. 


Hence,  we  conclude  from  (60)  that  cmtl/an_i,i  wil1  aPProach  zero  as  n*®’  as  long 
as  the  next  term,  <*,+itl/°*i,l»  also  approaches  zero.  Therefore,  the  approach  we 
shall  take  for  finding  the  is  to  assume  «n+1  .l/^N.l10  for  some  suitably 

large  N  and  then  use  (60)  to  compute  anj/otn-i,!  iteratively,  starting  at  n=N 
and  working  down.  This  procedure  yields  ai,i/ao,l  as  a  continued  fraction: 


«UL 

<*0,1 


YN 


(63) 


where  the  continued  fraction  terminates  at  yn* 
may  be  found  for  ^r0m  (61)*  However 

Y-n  =  Yn  « 


A  similar  continued  fraction 
,  since 

(64) 


the  result  is  just 


<*0,1  <*0,1 1 


(65) 


Hence  it  is  sufficient  to  evaluate  (63)  for  oq.,i/«o,l  and  then  use  (65)  to  find 


Suppose  the  continued  fraction  in  (63)  yields  a  value  r: 


We  may  then  substitute 

<*1,1  '■  ra0,l» 
and 

*-1,1  =  r**0,l» 


(66) 


(67) 

(68) 
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into  (51)  to  obtain 


[Bll+(l/4)£i  |Ui*|  {r+r*)]ao>i+ei2Qtn,2=*^del^lJd*‘ 

Solving  (69)  and  (52)  simultaneously  yields 

OQ  1  =  _ {Mi  &2?$\  turi*,- £l  2-e2.tud*,) _ . 

[Sn  +  (l/2)Ci  !Ui*j2Re(r)]622-|  B12I  2 

Once  oqj  has  been  computed,  oqj2  may  be  found  from  (52): 


<*0,2  = 


(69) 

(70) 


(71) 


We  have  now  developed  all  the  equations  needed  to  compute  the  periodic 
array  weights  in  Eq.(28).  To  summarize,  the  procedure  is  as  follows.  Starting 
from  the  signal  parameters  9^,  Cd»  9i »  £i »  <Mn'»  we  first  compute  the  vectors  Uq 
and  U-j  in  (7)  and  (14),  then  the  fljk  in  (42 )- (47 )  and  the  yn  in  (59).  From  the 
Yn,  we  evaluate  r  in  (66)  using  the  continued  fraction  in  (63).  (The  number  of 
terms  used  in  the  continued  fraction,  N,  is  chosen  large  enough  that  the  value 
of  r  obtained  is  insensitive  to  small  changes  in  N.)  After  r  is  found,  we  use 
(70)  to  find  aoj,  (71)  to  find  oq>2,  (67)  to  find  cqj,  (57)  to  find  ant\  for 
n*0,l,  and  finally  (58)  to  find  01^2  for  n*0.  From  the  an>i  and  «n>2  (and  using 
0^3=0),  the  Cn  may  be  evaluated  from  (36)  and  then  W(t')  from  (28). 

The  time-varying  weights  in  (28)  have  two  effects  on  array  performance. 
First,  they  cause  the  array  to  modulate  the  desired  signal.  (The  array  becomes 
a  time-varying,  or  frequency  dispersive,  channel [10]. )  Second,  the  array  output 
signal-to-interference-plus-noise  ratio  (SINR)  varies  periodically  with  time. 

Our  major  purpose  in  this  paper  is  to  evaluate  these  two  effects. 

Given  a  time-varying  weight  vector  W(t’),  the  desired  signal  component  of 
the  array  output  is 

Sd(t')  =  AdWT(t,)Udej(a«,t,  +  l(;d),  (72) 
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(73) 


(where  (iQ^^/kc2).  To  study  the  modulation  on  sd(t),  we  define 
ad(f)ejnd(t,)  =  AdwT(f)Ud. 

Then  ad(t' ) =Ad | ( t 1 )Ud|  is  the  envelope  modulation  and  nd(t 1 )=<W^(t ' )Ud  is  the 
phase  modulation.  Furthermore,  we  define  adn(t‘)  to  be  the  envelope  normalized 

to  its  value  in  the  absence  of  interference,  i.e.,  1 

i; 

adn(C)  -  ,•  <74>  ' 

Ad  I w0Tud | 

where  W0  is  the  steady-state  weight  vector  that  would  occur  without  1 

interference, 

W0  =  ($d  +  <t>n )  1S.  (75) 

(4>d,  On  and  S  are  given  in  (19), (21)  and  (23).)  We  present  our  results  in  terms 
of  adn(t')  rather  than  ad(t'),  because  the  effect  of  the  interference  can  be 
determined  most  easily  by  comparing  the  value  of  adn(t')  with  unity. 

The  output  desired  signal  power  is 

Pd(f)  =(  1/2)E  { |  sd(t 1 )  |  2  }  =  (l/2)Ad2|wT(t1)Ud|2.  (76) 

The  output  interference  signal  is 

Si(t')  =  Aicos(affl,f)ej(<^,t,  +  ','i)wT(t>)Ui,  (77) 

and  the  output  interference  power  is 

Pi(t')  =  (1/2)E  { |  s -j  ( t ' )  |  2 }  =  ( 1/2 )  A-j  2cos2(un)‘ t 1 )  |  W^(t 1  )Ud  |  2.  (78) 

The  output  thermal  noise  power  is 

2 

Pn(f)  =  f-  Wt(t')W(f).  (79) 

From  these  quantities,  the  output  interference-to-noise  ratio  (1NR), 
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and  the  output  signal-to-lnterference-plus-noise  ratio  (SINR), 

SINR  =  Pi(t^+7n(t')*  (81) 

may  be  computed. 

With  this  background,  we  now  discuss  the  numerical  results  obtained  from 
these  equations. 

III.  RESULTS 

In  this  section,  we  describe  the  effect  of  the  amplitude  modulated 
interference  signal  on  the  array.  In  Part  A,  we  show  typical  curves  of  desired 
signal  modulation,  output  INR  and  SINR  as  functions  of  time.  In  Parts  B-E,  we 
describe  the  effect  of  each  signal  parameter  on  the  desired  signal  modulation. 

In  Part  F,  we  assume  the  array  is  used  in  a  digital  communication  system  and 
show  how  the  bit  error  probability  is  affected  by  the  interference. 

A.  Typical  Waveforms 

Figure  3  shows  a  typical  curve  of  the  normalized  envelope  modulation 
adn^')  for  the  case  0d=O°,  0-j=5°,  Cd=10  dB,  =20  dB  and  fm'  =tHTi ' /2tt=2.*  As 
may  be  seen,  there  is  substantial  envelope  modulation  produced  by  the 
interference  for  this  set  of  parameters. 

When  the  phase  modulation  nd(t')  is  calculated,  one  discovers  an 

Interesting  result:  nd(t')  is  constant,  not  only  for  the  parameters  used  in 

Figure  3,  but  for  all  values  of  the  signal  parameters.  Thus,  this 

Interference  does  not  produce  phase  modulation  on  the  desired  signal.  A 

similar  result  was  also  found  previously  for  a  pulsed  interference  signal [4], 

With  pulsed  Interference,  It  Is  possible  to  prove  mathematically  that  there  is 

no  phase  modulation.  For  the  interference  considered  here,  we  have  not  been 

able  to  obtain  such  a  mathematical  proof.  However,  calculations  of  nq(t') 

based  on  the  equations  in  Section  II  consistently  show  that  there  is  no  phase 

*Because  the  frequency  of  the  periodic  term  in  (27)  is  2^',  the  envelope 
modulation  Is  periodic  with  period  l/2fm'=.25.  Figure  3  shows  one  period  of  the 
modulation. 
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modulation,  regardless  of  the  values  of  9^ ,  Ci  or  fm'.  This  result  is  of 
considerable  Importance  for  adaptive  arrays  to  be  used  in  communication 
systems  with  phase  modulation. 

Next,  Figures  4  and  5  show  the  output  INR  and  SINR  as  functions  of  time, 
over  one  period,  for  the  same  signal  parameters  as  in  Figure  3.  The  sharp 
drop  in  INR  in  Figure  4  at  t *  = . 125  occurs  when  the  incoming  interference  power 
is  zero.  The  variation  in  the  SINR  in  Figure  5  is  primarily  due  to  the 
variation  in  desired  signal  amplitude  as  shown  in  Figure  3.  However,  the  peak 
SINR  occurs  slightly  earlier  than  the  peak  of  adrift'),  reflecting  the 
influence  of  the  time-varying  interference  and  noise  powers  on  the  SINR. 

Figures  3,  4  and  5  are  intended  merely  to  illustrate  typical  array  behavior 
with  the  interference  modulation  in  Eq.(ll).  In  general,  the  desired  signal 
modulation  and  the  SINR  variation  change  substantially  as  the  signal  parameters 
0d,  £d»  0i »  Cl  and  fm‘  are  varied.  In  Sections  B-E,  we  describe  in  detail  the 
effect  of  each  signal  parameter  on  the  desired  signal  modulation.  In  Section  F, 
we  show  how  the  SINR  variations  affect  bit  error  probability  when  the  array  is 
used  in  a  digital  communication  system. 

To  describe  the  desired  signal  modulation,  we  shall  define  three 
quantities.  First,  we  let  amiq  and  amax  be  the  minimum  and  maximum  values  of 
adn(t')  during  the  modulation  period.  Then,  we  define 

m  =  >ax-.-a-rnln.  (82) 

amax 

m  is  the  envelope  variation  normalized  to  its  peak.  It  may  be  thought  of  as 
"fractional  modulation",  analogous  to  percentage  modulation  in  AM.  In  Parts 
B-E  below,  we  describe  the  effect  of  each  signal  parameter  on  amax  and  m. 

B.  The  Effect  of  Angle  of  Arrival 

Desired  signal  modulation  effects  are  small  unless  e-j  is  close  to  9<j. 

When  9i  is  far  from  9cj>  the  envelope  variation  m  is  small  and  the  peak  amax  iS 
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Figure  5.  S I NR  versus  time 

9d=rt°,  0i=5°,  5(1=1° 
Cl =20  dB,  fm'=2 
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close  to  unity.* 

Figures  6  and  7  show  typical  curves  of  m  and  a^y  as  functions  of  e-j  for 
9d=0°*  Sd=lf)  dB  and  Four  different  curves  are  shown,  for  r,i=5,  10,  IB 

and  20  d8.  It  is  seen  that  m  is  large  and  amax  is  small  only  when  9^  is  close 
to  0^.  Similar  curves  for  od*0  show  that  the  modulation  is  always  small 
unless  e-j^d. 

C .  The  Effect  of  Modulation  Frequency 

The  variation  m  and  the  peak  amax  are  large  at  low  fm‘  and  drop  as  fm‘ 
increases.  This  result  is  illustrated  in  Figures  8  and  9,  which  show  m  and 
amax  as  a  function  of  frequency  fm'  for  9d=0°,  £d=0  dB,  =20  dB,  and  for 
several  values  of  9-j . 

This  behavior  is  due  to  the  array  speed  of  response.  For  low  fm‘,  the 
array  is  fast  enough  to  follow  the  interference  modulation.  The  array  pattern 
changes  continuously  as  the  interference  power  varies.  At  high  f m * ,  however, 
the  array  weights  are  too  slow  to  keep  up  with  the  modulation.  Hence,  as  fm* 
is  increased,  the  array  weights  settle  into  steady-state  values  reflecting  the 
time-averaged  interference  power.  As  they  do  so,  the  desired  signal 
modulation  drops. 

Figure  8  shows  a  small  dip  in  m  for  intermediate  values  of  fm‘  for  some 

6-j .  This  effect  may  be  seen  in  the  6^  =  15°  curve,  for  example.  Such  behavior 

is  due  to  the  way  the  shape  of  the  modulation  envelope  changes  as  fm'  varies. 

Figure  10  illustrates  this  effect.  It  shows  a(jn(t')  as  a  function  of  time  for 

9j-15°  and  for  two  values  of  fm',  .015  and  .04.  (In  order  to  plot  the  two 

curves  on  the  same  graph,  we  have  plotted  both  versus  fm't'.)  Figure  10  shows 

that  a^x-aro-jn  is  essentially  the  same  for  f m  1  = . 04  as  for  fm'  =  .015.  However, 

m  is  lower  at  f ,n '  =  .015  than  at  f,n‘  =  .04  because  amax  is  higher. 

*An  exception  to  this  occurs  if  the  interference  causes  a  grating  null  in  the 
desired  signal  di rection[ll].  For  isotropic  elements  a  half  wavelength  apart  as 
considered  here,  there  are  no  grating  nulls. 
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6.  m  versus  e< 
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Figure  8.  m  versus  fm' 
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Figure  10.  a^d' )  versus  time 

«d*»Vei*15u,  TpT  dB,  ei  =20  dB 
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D.  The  Effect  of  Interference-to-Noi se  Ratio 

For  low  fm',  the  variation  m  is  largest  for  high  INR.  For  intermediate 
values  of  fm',  m  peaks  at  intermediate  INR.  amax  is  unity  at  low  fm‘  and  drops 
to  a  minimum  value  as  fm'  increases.  The  higher  the  INR,  the  farther  amax  drops 
at  high  fm‘.  Figures  11  and  12  illustrate  these  results  for  the  case  0^=0° , 
0i=5°,  £^=0  dB  and  £i  between  0  and  30  dB. 

For  intermediate  fm‘,  it  may  be  seen  in  Figure  11  that  m  peaks  at 
intermediate  INR.  This  behavior  is  due  to  the  way  the  modulation  envelope 
changes  with  INR  and  to  the  way  m  is  defined.  To  illustrate  this.  Figure  13 
shows  the  modulation  envelope  a(jn(t')  over  one  period  for  f m ' =1  and  several 
INR.  It  may  be  seen  that  amax-am-jn  at  first  increases  with  INR  but  then 
decreases,  because  amin,  the  smallest  value  of  afjn(t')  during  the  cycle,  does 
not  drop  below  about  0.17. 

E.  The  Effect  of  Desi red  Signal-to-Noise  Ratio 

The  variation  m  is  largest  and  the  peak  amax  is  smallest  for  low  r,^.  As 
is  increased,  m  decreases  and  amax  increases. 

Figures  14  and  15  illustrate  this  behavior.  Figure  14  shows  m  and  Figure 
15  shows  amax,  both  versus  fm'  for  9^=0°,  6i=5°,  Ci =20  dB  and  for  five  values 
of  £d  between  0  and  20  dB.  It  is  seen  that  for  a  given  SNR,  m  peaks  slightly 
at  intermediate  values  of  fm'.  This  behavior  is  due  to  the  changes  in  the 
shape  of  the  modulation  envelope.  As  an  example.  Figure  16  shows  a(jn(t') 
versus  t'  for  r>C|=0  dB  and  10  dB  and  for  fin’  =  .01,.12  and  10.  These  values  of 
fm'  span  the  peak  of  the  curve  for  f>(j=10  dB  in  Figure  14  and  show  how  the 
shape  of  ( t  ’ )  changes  with  fm'.  For  r,d=0  dB,  the  curve  of  m  vs.  fm'  in 
Figure  14  shows  no  peak,  and  this  behavior  is  seen  in  Figure  16. 

F .  Bit  Error  Probabi 1 ity 

The  effect  of  a  time-varying  SINR  can  be  evaluated  most  meaningfully  by 

computing  the  bit  error  probability  when  the  array  is  used  in  a  digital 
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Figure  11.  m  versus  fm* 

fld=7ia;“ei'=5n,  ?d=0  dB 
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communication  system.  For  this  purpose  we  arbitrarily  assume  the  desired  signal 
to  be  a  differential  phase-shift  keyed  (DPSK)  biphase  modulated  signal[12]. 


Also,  we  assume  the  reference  signal  is  a  replica  of  the  desired  signal.* 

Because  the  interference  does  not  produce  phase  modulation  on  the  desired 
signal,  the  only  effect  of  the  interference  on  system  performance  will  be  to 
vary  the  array  output  S I NR  and  hence  change  the  bit  error  probability.  For 
a  DPSK  signal  in  white  noise,  the  bit  error  probability  Pe  is  given  by [14]: 

Pe  =  (l/2)e'Eb/No,  (83) 

where  Ep  is  the  signal  energy  per  bit  and  N0  is  the  one-sided  thermal  noise 
spectral  density.  For  our  purposes,  it  is  convenient  to  express  Eb/N0  in  terms 
of  signal -to-noise  ratio.  Since  Ep=PdTp,  where  P<j  is  signal  power  and  Tp  is  the 
bit  duration,  and  since  1/Tp  is  the  effective  noise  bandwidth,  N0/Tp  is  the 
received  noise  power.  Hence, 

*0  ■  tw  ■ SNR- 

In  addition,  for  this  analysis  we  shall  assume  the  interference  power  at  the 
array  output  has  the  same  effect  on  detector  performance  as  thermal  noise 
power.  With  this  assumption,  the  bit  error  probability  may  be  written 

Pe  =  (l/2)e'SINR.  (85) 

Finally,  since  the  interference  causes  the  SINR  to  vary  periodically,  we  obtain 
the  effective  bit  error  probability  by  averaging  Pe  over  one  period  of  the 
interference  modulation: 


7*  *  2fn'/<1/2fm')(l/2)e-SINR(t')df. 
0 


(86) 


*As  long  as  the  desired  signal  bandwidth  is  not  excessive,  these  assumptions 
produce  the  same  ^  as  in  (19)  and  the  same  S  as  in  (23),  and  hence  will  yield 
the  same  weight  behavior  as  the  CW  signal  assumed  above.  Moreover,  it  has  been 
shown  that  desired  signal  bandwidth  has  almost  no  effect  on  array  performance 
even  if  the  bandwidth  is  large[13]. 


This  calculation  is  valid  as  long  as  the  bit  rate  is  large  compared  with  the 
interference  modulation  period,  which  we  assume  to  be  the  case. 

Figure  17  shows  typical  curves  of  "Pe  as  a  function  of  fm'  for  9^=0°, 

0i=3O°,  ^=6  dB  and  for  several  values  of  between  -20  dB  and  30  dB.  Figure 
18  shows  "Pe  versus  fm‘  for  9(1=0°,  ^=6  dB,  r,-j=20  dB  and  for  several  9i  between 
5°  and  20°.  Both  these  figures  show  that  fm'  has  almost  no  effect  on  bit  error 
probability.  At  low  fm‘,  where  the  weights  can  track  the  interference 
modulation,  the  SINR  varies  over  a  wider  range  than  at  high  fm',  where  the 
weights  are  too  slow  to  follow  the  modulation.  Nevertheless,  the  average  bit 
error  probability  is  essentially  the  same  for  all  fm'. 

The  curves  in  Figures  17  and  18  show  a  slight  increase  in  as  fm‘ 
increases.  At  high  fm‘,  where  the  array  is  too  slow  to  track  the  interference 
power  variation,  the  weights  are  constant.  Their  values  are  determined  by  the 
time-average  interference  power,  Aj2/2  (see  Eg. (20)),  or  equivalently  by  the 
time-average  input  INR,  Zj/2.  In  this  case  one  can  show  that  Pg  has  the  same 
value  as  the  bit  error  probability  for  an  array  receiving  CW  interference  with 
input  INR  of  K\/2.  ~Fe  is  slightly  lower  at  low  fin'  because  in  this  case  the 
weights  track  the  modulation  and  yield  maximum  available  SINR  for  each 
instantaneous  value  of  interference  power.  However,  the  difference  between  "Pe 
at  low  and  high  fm‘  is  so  small  that,  for  practical  purposes,  "Pg  may  be 
considered  constant,  equal  to  that  for  CW  interference  with  INR  =  r,-j/2. 

IV.  CONCLUSIONS 

We  have  examined  the  effect  of  a  sinusoidally  modulated  (double-sideband, 
suppressed-carrier)  interference  signal  on  the  performance  of  an  adaptive  array. 
The  major  effect  of  such  interference  is  to  cause  envelope  modulation  on  the 
desired  signal.  This  envelope  modulation  is  large  only  when  the  interference 
arrival  angle  is  close  to  that  of  the  desired  signal  and  when  the  interference 
modulation  frequency  is  low  relative  to  the  array  speed  of  response.  The 


Figure  17.  Bit  Error  Probabi 1 i ty  versus^  fml 
■9^0^97=30°, 
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envelope  modulation  is  greatest  when  the  interference  power  is  high  and  the 
desired  signal  power  is  low. 

The  modulated  interference  produces  no  phase  modulation  on  the  desired 
signal.  Moreover,  if  the  desired  signal  is  a  digital  communication  signal,  such 
interference  results  in  a  bit  error  probability  essentially  the  same  as  that  for 
CW  interference  (from  the  same  angle  and  with  power  equal  to  the  time-average 
power  of  the  modulated  interference).  Interference  modulation  frequency  has 
almost  no  effect  on  bit  error  probability. 
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